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SUMMARY 


A transformation theory of systems of partial differential 
equations is developed which allows the construction of classes of 
pressure— density relations depending on arbitrarily many parameters 
for which the equations governing the flow can be transformed into an 
essentially simpler form., namely, into the Cauchy— Riemann equations in 
the subsonic region, into the system corresponding to the wave equation 
in the supersonic case, and finally into that form corresponding to the 
Tricomi equation in the transonic region. The transition from one 
system of differential equations to the other is always such that only 
the solving of ordinary differential equations is required in order to 
find solutions of the more complicated system from corresponding 
solutions of the simpler one. 


INTRODUCTION 


It has already been observed by Chaplygin (reference 1 ) that in 
the case of a pressure— density relation of the form p = — — + b,a > 0 

the equations describing a steady irrotational flow in the plane can be 
reduced to the Cauchy-Riemann differential equations. This is achieved 
by transforming the equations into the hodograph plane, that is, the 
plane of the velocity components , and using suitable combinations of 
the latter as independent variables. Since any pressure-density 
relation can, in smaller regions of the variables involved, be approxi- 
mated by the foregoing relation, approximation theories can be developed 
which make 1 use of the Cauchy-Riemann equations. A series of recent 
papers, by Von Harman, Tsien, Bers, Gelbart, Bartnoff, and Lin, follow 
this idea (references 2 to 9)» 

The investigations taken up in this report originated in the 
question whether the ideas developed by the Swedish mathematician 
Baecklund, in his transformation theory of partial differential equa- 
tions of second order, could, after suitable modifications, be utilized 
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for the study of the equations of compressible fluid flow. The funda- 
mental problem of Baecklund was to find all pairs of surfaces Z 
and Z* in an x,y,z— space and an X* ,y* , z’— space, respectively, and a 
one—to—one mapping between them such that, for corresponding points, 
four given equations 


^fcy^jPj^xSySzSpSg.') =0 (i = 1,2,3, A) (1) 


are satisfied. As usual p and q designate the partial derivatives 
of z with respect to x and y and, similarly, p’ and q* the 
derivatives of z* with respect to x’ and y*. It is clear that, 
for a surface Z given, in general no surface Z * can be found which 
can be mapped onto Z in such a way that equations (1) are satisfied. 
In order to find the necessary and sufficient conditions for the proper 
surface I* and mapping of Z onto Z ! to exist, the following 
operations, according to Baecklund, must be performed (reference 10): 

(a) Considering x* and y* as functions of x and y, each of 
equations (1) is to be differentiated up to the second order with 
respect to the two independent variables. 

(b) liom the equations thus obtained combined with equations (1), 
all quantities bearing a prime are to be eliminated. 

In general the elimination result will consist of two partial 
differential equations of the third order in z(x,y). They represent 
the conditions on the surface Z . If the differentiations of first 
order in step (a) combined with equations (l) already allow for the 
elimination of the primed quantities, the resulting equation of second ' 
order represents the conditions on the surface Z sought. Steps (a) 
and (b) may also be taken with reversed roles of surfaces Z and Z* 
and will lead to one or two partial differential equations in z* as 
a function of x* and y* . 

These considerations show that equations (l) link solutions of 
two, in general, different systems of partial differential equations 
with each other. If one system is simpler than the other, assertions 
about the solutions of the more difficult system may be derived from 
knowledge of the solutions of the simpler one. 

The steady irrotational flow of a compressible fluid may be 
described by the equations 

Tiy - Tx = 0 


(pu) x + (pv) y = 0 


( 2 ) 
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where x and y are Cartesian coordinates in the physical plane 
and u and v the components of the velocity vector. The density p 
is to "be th ought of , according to Bernoulli's lav, as a 'function of the 

speed q = \/u 2 + v2. The subscripts x and y in system (2) designate 
the respective differentiations. 

Mother description of the flow is obtained by introducing the 
potential function f and the stream function as follows; 

Equations (2), considered as integrability conditions, show that there 
exist functions 0(x,y) and i|r(x,y) such that 

u = 0x P u = Vy 

v = fy -pv = tjj. 

This leads to the equations 


*X = “P0y 


*y » P0 X 


(3) 


in which the arguments u and v in p are to be replaced by f x 
and <f>y, respectively. 

Both systems (2) and (3) are nonlinear in the case of compressible 
fluid flow because p is a variable function of u and v. As was 
recognized by Chaplygin, systems (2) and (3) take on essentially simpler 
form if the physical plane is replaced by the hodograph plane; that is, 
if the velocity components u and v are introduced as independent 
variables instead of x and y. Systems (2) and (3) then become linear. 
Each of the systems takes on the form 


Fi I. 


u 


+ °i 5, 


4 


T). + 

U 


4 


T) = 0 (i = 1,2) 




where g and q represent the coordinates x and y in the case of 
equations (2) and the potential and stream functions } and i{r in the 
case of equations (3). 



4 


MCA TI'T 2065 


1 


The coefficients F^. and are definite functions of u and v 

and their forms depend on the equation of state adopted. For an incom- 
pressible fluid, equations (4) are identical with the Cauchy-Bi emann 
equations . 


Very often it is convenient to introduce coordinates other than u 
and v in the holograph plane, for example, polar coordinates, or any 
other generalized coordinates s and t connected with u and v 
hy a transformation 


or 


s = s(u,v) 
t = t(u,v) 


(5) 


u = u(s,t) 
V = v(s,t) 


(5a) 


where "both transformation (5) and its inverse, transformation (5a), are 
continuously differentiable. If such a transformation is applied to 
equations (4), a system 


f l 5 + Si §. + fo 11 + 4 \ = 0 C 1 “ 

x s x t ^ S ^ t- 


(6) 


of the same form as system (4) is obtained with new coefficients f^ 
and gj£ which are now to "be written as functions of the new coordi- 
nates s and t. 

Following the ideas of Baecklund, the following question may be 
raised: Is it possible to interrelate two given systems of differential 

equations of the form of system (6) by a transformation of the Baecklund 
type? As an analogue to equations (1), a transformation of the 
“Baec kl und type" will be defined here as given by a system of equations 
of the form 

^i( s ,t,| ,T),| s ,|t ■ ,T ls jT it' ,s * * *5 * s 1 * t * ,r ^ *s ’ *t 0 = ^ ^7) 


(i = 1,2, . . .,6) 
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A solution of system (7) is represented by two pairs of functions 

i = 6(s,t) 

T\ = Tj(Sj t ) 

having as independent variables s,t and s^t', respectively, and a 
one-to-one mapping between the s ,t— and s’, t’— plane such that at corre- 
sponding points all equations (7) are satisfied. 

If a geometrical interpretation is desired, two spaces P and P’, 
each of four dimensions, with the coordinates (s,t,g,q) 
and (s’,t’, l*,!}* ), respectively, have to be introduced. A pair of 
functions g = g(s,t) and T] = q(s,t) represents a two-dimensional 
surface Z in P and correspondingly a pair of functions £* = g* (s*,t* ) 
and tj* = T)*(s*,t*) represents a surface 2* in p*. To solve 
system (7) means to find all pairs of surfaces 2 and Z* and a 
one-to-one mapping between them such that the quantities character- 
izing their tangent planes at corresponding points satisfy 
relations (7). 

That the choice of six equations in system (7) is natural follows 
from the fact that if s and t are considered, as independent vari- 
ables and £,ti,£ , , t 1 , .> s, .> t* as dependent variables, the number of 
equations coincides with the number of unknown functions. 

Again, as in the case of equations (1), the following important 
question presents itself: What surfaces Z (or 2*) belong to a solu- 

tion of system (7)? It is to be eapected that necessary and sufficient 
conditions may be found from operations corresponding to those applied 
to equations (l), that is, 

(a) Considering s’ and t* as functions of s and t, each of 
equations (7) is to be differentiated with respect to s and t Tip to 
the second order. 

(b) All quantities bearing a prime are to be eliminated from the 
equations thus obtained combined with equations (7). 

The. result of the elimination will in general consist of four 
partial differential equations of third order. In particular cases a 
differentiation of first order in operation (a) suffices to yield 
necessary and sufficient conditions on Z so that Z’ should exist, 
by the elimination of the primed quantities. This is the case with 
all Baecklund transformations studied and applied in this report. 


I* = 6*(s*.,t») 

V = T)‘(s*,t*) 


(8) 
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The essential point in the preceding consideration is again the 
fact that, by a Baecklund transformation of the type of system ( 7 ) , two, 
in general, different systems of partial differential equations con- 
taining two independent and two dependent Yariahles are interconnected. 


For the purposes of this investigation it is not necessary to 
study the most general type of Baecklund transformation such as 
system ( 7 ) . Since the interest is directed to linear systems of the 
type of system (6), it is sufficient to consider only such equations ( 7 ) 
which are linear and homogeneous in the quantities and q * and 

their derivatives of first order. The further more restrictive assump- 
tion will "be made that two of equations (7) describe a one-to-one 
mapping of a domain of the. s,t— plane onto a domain of the s‘,t’— plane 
or, in other words, represent a transformation of the independent vari- 
ables given by 


s* = f(s,t) 

t* = g(Sjt) 


(9) 


Since transformations of any quantities resulting from such a 
coordinate transformation can easily be handled, there is no essential 
loss of generality if the assumption is henceforth made that equa- 
tions (9) reduce to 



( 10 ) 


System (7) can then be replaced by a system of only four equations, 
linear and homogeneous in | ,q,£ ’ , and q*, and their derivatives of 
first order with respect to s and t. The coefficients may be arbi- 
trary functions of s and t satisfying the condition that they be 
twice continuously differentiable. Only the additional condition, that 
the four equations can be solved for the four first partial derivatives 
of | and q and also for the corresponding derivatives of § * 
and q’, is assumed. Equations (7) can therefore be brought into the 
form 
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T) + a,l s 
‘t 1 & 


+ a £ 


n + 


i’ 


+ y\ V 


*'t = CT i* 


+ cr^" T1 + CO"*" t 

s 2 's 1 § 


+ CX?~ T1 + 

t 2 't 


p i 5 + 




+ 6 1»’ 


R 1 » 
6 2 11 



PT S, 



TJ + S + T1 + Ct£ g 
S 1 5 t 2 't 1 5 


+ a 2 V+^' 


+ y\ *i* 


(11) 


i*t 


= a i is + 4 % + “ 1 6 t + 4 \ + p i 6 + 4 * + 4 i‘ + 4 v 


where the coefficients on the right are all functions of s and 
and the determinant of the coefficient matrix of equations (11) 
contains d in the first four columns on the right is supposed to "be 
different from zero. 


Equations (11) can he brought into a sinpler form by use of the 
matrix calculus notation. Set 
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Further designate the column matrix by £ and (^ 1 ) "by £». if 

differentiation of a matrix means, as usual, differentiation of each of 
its elements, equations ( 11 ) can then he written in the simple form 


^*S = Hs + nt + H + C Z' 
t't = S ^s +B ^ + - 


( 12 ) 


A systematic study will now he made of the possibilities of transforming 
systems of equations, of the form of system ( 6 ) into each other hy trans- 
formations of the type described hy equations ( 12 ), and applications to 
compressible fluid flow will he made.- 

This investigation was conducted at Syracuse University under the 
sponsorship and with the financial assistance of the National Advisory 
Committee for Aeronautics. 

The author is indebted to Messrs. Helmut Aulbach and Samuel 
Schechter for valuable assistance. 


SYMBOLS 


a,b,c 


i R i i R i 


pScr 1 /!- 3 -,^ 3 - 
k k k k 


constants 


functions of s and t defined by equations (ll) 


A,B,C,D,K,S,T,W matrices having the previous functions as respective 

elements 

F-j^G^ functions defined by equations (l) and ( 7 ), 

respectively 


H 


matrix defined by equation (lh) 


I 


unit matrix 
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K 

^( 8 )^( 8 ) 

k(ff) 

M 

P 

q. 

s,t,s*,t* 

u,v 

w - 1 

¥* 


ajP 


7 

4 

S^pW ^ 1 

0 ■ 


i 

p 


a 


T 

9 

K 

t,t* 

cn 


constant matrix 

functions defined in theorem. 3 (appendix) 

major ant function used in theorem 1 (appendix) 

Mach number 

pressure 

speed 

curvilinear coordinates in holograph plane 
velocity components 

Cartesian coordinates in physical plane 

inverse of matrix W 
"adjoint" of matrix V 
constants 

exponent in pressure-density relation 
elements of the matrix H 


functions defined by equations ( 35 a) 
velocity potential 
stream function 
density 

function defined by equation ( 87 ) 
variable defined by equation ( 80 ) 
angle between velocity vector and x-axis 
auxiliary variable defined by relation ( 74 ) 
column vectors 




K* 1 / W/ respectivel3r 

function defined by equations (49) 
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o>*(s*) function defined "by formula (6l) 

t(m) trace of matrix M 

functions of the form given "by equation ( ill ) which 
J approximate ou 

cu*^_,(n* 2 functions which approximate cd* 

s as subscript, partial derivative with respect to s; 

any other variable used as subscript, partial 
derivative with respect to that variable 


ANALYSE 

1. Baecklund Transformations Connecting Two Systems of 
Partial- Differential Equations of First Order 


A steady irrotational flow in the plane is described by a system 
of differential equations of the form of system (6). A study of the 
application of Baecklund transformations to such systems must start 
with answering the following fundamental question: Which of the 

Baecklund transformations of the form of equations (12) transform two 
systems of partial differential equations of the form of system (6) 
into each other? A complete answer will be given in this section, and 
it will appear in the form of a system of partial differential equa- 
tions in the matrix coefficients in equations (12). 

In order to simplify the writing, system (6) will always be 
written in a matrix form. After setting 



system (6) may be written in the form 


+ G£t = 0 . 


( 13 ) 
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(Henceforth, a system of equations -written in matr ix form as a s ing le 
relation will also "be referred to as an equation. ) Without loss of 
generality it may "be assumed that the determinant of F is never zero. 
As can he easily seen, this means that a line t = Constant can never 
he touched hy a characteristic curve of equation ( 13 ). This can always 
he achieved hy a suitable choice of the coordinates s and t. 

Equation (13) can now he written in the form 


£s=H£ t (1M 

where H = — E - ^G- represents a matrix function of s and t. 

In order to find the conditions £ must satisfy if it is to 
belong to a solution of system (12), the latter equations must he 
differentiated, and the primed quantities eliminated. If the first 
equation is differentiated with respect to t and the second with 
respect to s, subtraction of the ensuing equations yields 

0 = S? ss + (W — E)£ st - T5 tt + 


(S s - B t + B)5 S + (W s - T t - A) ?t + (B s - + 

B£*« " Cft + ( B s - °t)V ( 15 ) 

If, further, the expressions for £ 1 s and £ \ from equations ( 12 ) 
are substituted in equation (15), the resulting equation is 

0 =S£ ss + (W-KK st -T£ tt + 

(S s - R t + B + DR - CS)^ S + (V s - T t - A + DT • 

^DA - CB + Bs - A t )g + (DC - CD + D s - C t ^ * 

If the assumption is introduced that the coefficient of 
tion (l6) is identically zero, or 

D s - C t + DC - CD = 0 (17) 

then equation (16) is already free of primed quantities and therefore 
represents a condition on {; . In general, equation (16) is of second 


- cw)e t + 

(16) 

£ 1 in equa— 
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order. Since only equations of the form of equation (13) are studied 
here, the further assumptions that 


W = B 
S = T = 0 


1 


(18) 


Bg — A^ + BA — GB — 0 

•will "be made. Transformation (12) then reduces to 

5*8 = vts + K + eg' 

g « t = Wgt + Bg + Dg» 


( 19 ) 


( 20 ) 


and equation (l6) to 

(-*t + B + DW)g s + (W s - A - CW)g t = 0 (21) 

The latter procedure leading to equation (21) may he subordinated 
to the derivation of the integrability conditions of a "Mayer— Lie 
system" of partial differential equations. (See reference 11. ) Indeed, 
if g is considered as given, then equations (20) represent such a 
system for the determination of g*. A Mayer— Lie system expresses all 
derivative's of first order of the unknown functions in terms of the 
independent and the dependent variables themselves. The principal 
theorem of the theory of Mayer— Lie systems applied here asserts that, 
if equations (17) and (19) are satisfied, equation (21) is not only a 
necessary but also a sufficient condition for the existence of a g * 
that can be connected with g by equations (20). It furthermore 
asserts that there exists a g* having an arbitrarily preassigned 
initial value g* 0 at a given point' (s 0 ,t 0 ), and g* is uniquely 

determined by g : q. The computation of g* requires the solving of a 
system of ordinary differential equations with given initial conditions. 

It is now easy to find further conditions on the coefficient of 
equations (12) so that also g* will satisfy a system of the form of 
equation (13). The roles of the quantities g and g* have to be 
reversed. To attain this, equations (20) have first to be solved 
for g s and g-j-. The assumption that this is possible was introduced 

from the begi nni ng and means, evidently, that the determinant of ¥ is 
nowhere zero or that \f has an inverse, W* = W - 1. The inversion 
yields 

g s = W*g' s + A’g* + C’g 1 


£t = w ’£’t + + 


( 22 ) 
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with 


l 


A* = -T^C, C* = -W^A 
B* = -W“^D, D» = -W^B 


W* = Tf 1 


(23) 


The conditions sought will now he obtained, by replacing all quantities 
in equations (17) and (19) by the corresponding primed quantities. 

This leads to 


D* s - C«t + D’C’ - C*D* = 0 (24) 

B« s - A‘ t + D*A» - C’B * = 0 (25) 


and the equation satisfied by {;* is, corresponding to equation (21), 
(-W‘ t + B* + D*W*)£‘ s + ^W’ s - A* - C»W*)g » t = 0 (26) 


The final result obtained may be stated as follows: If all the differ- 

ential equations (17), (19), (24), and (25) in A, B, C, D, and W are 
satisfied, the Baecklund transformation, equations (20), transforms 
system (21) into system (26). 

, If the system of equations consisting of equations (17), (19), 

(24), and (25) is written in scalar form, the number of unknown scalar 
functions exceeds the number of equations by four. It is therefore to 
be expected that, if four of the unknowns are arbitrarily chosen or, 
more generally, four additional relations between the u nkn owns are 
added to the system, there will still exist infinitely many solutions. 

In particular, it is to be expected that the equation in £ or £ * 
may arbitrarily be prescribed and still be transformable into infini tely 
many other equations by suitable Baecklund transformations. 

Consider the problem of finding all Baecklund transformations such 
that a preassigned equation 


5*8 = n’t 


(27) 
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implies equation (26). (Naturally, at all points (s,t ) where the 
coefficients of equation (26) are not hoth singular, that is, do not 
have vanishing determinants, equations (26) and (27) express the same 
relation.) A comparison of equations (27) and (26) shows that then 

W* s - A* - C*W* = — B* — D*¥*)h (28) 


This equation must he added to equations (17), (19 ) 3 (24), and (25). 
The whole system of equations can he essentially simplified hy the 
following steps: If the second of equations (20) is multiplied hy H 

from the left and then subtracted from the first, the relation 


0 = W£ s - EW£ t + (A - HB)£ + (C - HE>)£ * 


(29) 


is obtained. If the coefficients of equation (26) are not hoth singular, 
equation (29) will hold for any £ and £ * connected hy the Baecklund 
transformation, equations (20). An arbitrary constant can he chosen 
for £, and for a given point (s,t) a corresponding value of £ * can 
he preassigned. A substitution into equation (29) then shows that 


A = BB and 0 = HD 


(30) 


Tn the degenerate case where the coefficients of equation (26) are 
singular, equations (30) will he assumed to he true too. System (29) 
becomes now 


5 b 


(31) 


Thus, relations (30) and (28) are obtained and are to he added to 
equations (17) , (19 ) 3 (24), and (25). 

It will now he shown that equations (24) and (25) are already 
consequences of equations (17), (19) j (28), and (30). Indeed, equa- 
tions (23) imply that 


D« s = -N‘ S B - W*B S , C' t = -N« t A - W'A,. 
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A substitution into the left side L of equation (24) leads to 
L = D* S - C ! t + D»C* - C*I>* = -¥*gB + W ! t A - ¥ ! (B S - + ¥‘B¥ ! A - ¥‘A¥*B 

and by use of relation (19) 


L = — ¥*gB + ¥‘ t A + ¥* (DA - CB) + ¥‘B¥*A - ¥'.A¥ , B 
= ^¥« t - B* - D*¥*)a - (¥‘ s - A* - C*¥’)b 


But this expression is zero on account of equation (28) and the first 
of equations (30)- In a similar way it can be verified that equa- 
tion (25) is also a consequence of equations (17) j (19), (28) * and (30). 

It is easy to eliminate the primed quantities from equation (28) 
by direct substitution or., more simply, by observing that the new equa- 
tion may be obtained from equation (28) by dropping the primes from the 
latter and replacing H by ¥“lH¥ which corresponds to an interchange 
of the roles of £ and (; * . This leads to 

¥ s — A - C¥ = (¥ t - B - D¥^¥-lH¥ (28a) 


Equations (30) allow the elimination of A and C from equations (17), 
(19), and (28a). If this is done the following final system for B, D, 
and ¥ is obtained: 


D_ - (ED), + DHD - HD 2 = 0 

S u 

Bg — (HB)^. + DHB — HDB — 0 

¥ s - HB - HD¥ « (V t - B - D¥)w _1 E¥ 


(32) 


If equations (32) are satisfied, the Baecklund transformation described by 
equations (20), with A and C defined by equation (30), transforms the 
equation £ s = ¥^-E¥^ into the equation £* g = 
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S.ome general remarks 'will now "be made about tide nature of equa- 
tions (32)j The first equation contains only the unknown D, the 
second contains the unknown B in addition to D, and finally the 
third contains all three unknowns D, B, and ¥. If the equations 
are solved successively, one of the unknown matrices has to he deter- 
mined at every step. Since the whole system of equations (32) is of 
the Cauchy-Kowalewski type and of first order, the general solution 
will depend on 12 arbitrary functions of 1 variable which arise from 
the initial conditions. This shows that, in general, each of a mani- 
fold of systems of equations depending on 12 arbitrary functions of 
1 variable can be transformed into a given system. 

A second remark refers to the possibility of generalizing the 
transformations thus far considered. Transformations (20) do not form 
a group. This means that, if two such transformations are composed, 
the resulting transformation will in general not be contained among 
transformations (20), but it is exactly this fact that allows for the 
construction of an extended class of . transformations which make possible 
the linking ’ of systems of differential equations not transformable into 
each other before. By repeating the composition process several times, 
the class of transformations can be widened more and more and an 
increasing flexibility can then be achieved. 

The differential equations connected with the study of a steady 
irrotational flow are of a special type. As a consequence, it is not 
necessary to consider here the’ most general transformation of the form 
of equations (20). The proper restrictions to be made on these trans- 
formations will be discussed in the following section. 


2. Specialized Systems of Differential Equations and 
Specialized Baecklund Transformations 

The systems of differential equations £ s = occurring in the 

study of a steady irrotational flow have, after a suitable choice of 
coordinates s and t in the hodograph plane, the following two 
properties: (a) The elements in the principal diagonal of the 

matrix H are both identically zero and (b) H depends only on one of 
the variables s and t, say, s. 

Property (a) expresses the fact that the equations can be derived 
from a problem in the calculus of variations. 

Properties (a) and (b) suggest the consideration of only such of 
transformations (20) which preserve both properties (a) and (b). 
Property (b) is certa inly not destroyed if the assumption is made that 
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all coefficients Ixl transformation (20) depend only on s. This leads 
immediately to a remarkable simplification of system (32). The 
res alting system of ordinary differential equations is 


D s + DSD - HD 2 = 0 
B s + DHB - HDB = 0 


(33) 


W s — HB — HDW = — (B + DW )'W~ ^HW J 

The subscript s indicates the ordinary derivative with respect to the 
only independent variable s. 


Before conditions are introduced insuring preservation also of 
property (a), some consequences will be drawn from system (33). 


The following notation of matrix theory will frequently be used: 

mA 


/ m^ 

(1) The trace t(M) of a matrix M = j x 




4 ) 


is defined as 


the sum of the elements lying in its principal diagonal, that is. 


t(M) = m.1 + m 2 

The trace t is an additive function of M, that is, for any two 
matrices Mq and Mg, 

t (Mq + M 2 ) = t (Mq) + t (M 2 ) 
and the trace of a scalar multiple A,M of M is 

r(XM) = Xt (M) 


It is further easily verified that for a product of the matrices 
and Mg 


t (MqM2_) = T (M 2 Mq) 
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or that 


T (M]_M2 - M 2 M x ) = 0 


(2) The "adjoint" M* of a matrix M is defined as the matrix 



Evidently MM* = M*M = 8 (M)l, where I 


represents the unit matrix. 



and 


5(M) 



mP- 

2 

4 


is the determinant of M. 


If BOO * °, J£L = ri. 


The simple facts of matrix theory reviewed under definitions (l) 
and (2) will now he used to draw consequences from equations (33)- If 
the trace is taken of the left side of the first of equations (33), the 
following equation is obtained; 


t(D s ) = t(HD 2 - DHD) = t(HD'D - D-HD) = 0 


or 

t(D s ; = [t(D)] s = 0 


This shows that the trace of D is constant. 
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If the same equation is 
and the trace is then taken , 


multiplied, from the right side, by 
the resulting equation is 


D* 


t.(D s D*) = 5(D)t(HD - DH) = 0 


HoweYer, an easy confutation shows that 


t(D s D*) = ( 8 ^ sf - ( S l) s af - ( 8 |) s 5^ + (»|) o - [ 8 (B)], 

This leads to the fact that the determinant of D is also const an t. 


In a similar way, multiplying the second of equations ( 33 ) hy B* 
and the third hy W* from the right and computing the trace afterward 
lead, as the reader mil easily verify, to the further result that 
the determinants of B and W are constant too. 

If a solution of system ( 33 ) is known, a four— parameter family of 
new solutions may he obtained hy leaving D unchanged hut replacing B 
and W hy new matrices B = BK and W = WK with a common constant 
right-hand matrix_f actor _K, whose determinant 5(K) / 0. Indeed, if 
substitution of B and W in the second and third of equations ( 33 ) 
is made, the matrix K appears in each term as a common right-hand 
factor, and the equations obtained are therefore consequences of 
equations ( 33 )- 

A final remark concerns some particular solutions of equatioxis (33): 
If D represents a solution of the first of equations ( 33 ) 
and 8 (D) / 0, then any scalar multiples of D 


B = ID and W = pD 


with constant factors X and p, p / 0 , satisfy the second and third 
of equations (33)- The proof is immediately obtained hy substitution. 

The general solution of system ( 33 ) depends on 12 arbitrary 
constants, hut only those solutions which guarantee the preservation 
of property (a) are of interest. By the addition of this new postulate, 
a subclass of the solutions of system (33) will he obtained. Without 
trying to characterize the latter completely, a six— parameter family 
of solutions of this subclass will be constructed by a specialization 
of the form of the matrices involved. 
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A simple method of assuring the preservation of property (a) is 
to assume ¥ to he of the form 


¥ = 


(34) 


Indeed, the matrix H of the primed equation £* = HfV is supposed 

to have the form t 


H = 


and the matrix ¥“^B¥ of the transformed equation = ¥~ is 
then / 


¥~ 1 Bf = 


'il 

ox. 



which shows that property (a) is preserved. 

Tn order to obtain solutions of system ( 33 ) with the matrix ¥ ' of 
the form of equation (34) j an adjustment of the forms of T> and B 
must he made hy setting 


D = 


and 


B = 


0 


Pi 


4 


0 


(35) 


It will now he proved that system (33) has solutions of the form of 
equations ( 34 ) and ( 35 ) with arbitrarily preassigned initial values of 

12 12 1 2 

the six available matrix elements 5, 5, 8 , 8 , 03 , and co . 

2 1 2 * 1 J 1 J 2 
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If the same equation is multiplied, from the right side, by D* 
and the trace is then taken, the resulting equation is 


t(D s D*) = 5(D)t(HD - DH) = 0 


However, an easy computation shows that 


T(D S B*) = ( 5 ^ 




Sj = [8(D)] S 


This leads to the fact that the determinant of D is also constant. 


In a similar way, multiplying the second of equations ( 33 ) by B* 
and the third by ¥* from the right and computing the trace afterward 
lead, as the reader will easily verify, to the further result that ■ 
the determinants of B and ¥ are constant too. 


If a solution of system (33) i s known, a four— parameter family of 
new solutions may be obtained by leaving D unchanged but replacing B 
and ¥ by new matrices B = BK and ¥ = WK with a common constant 
right-hand matrix_f actor _K, whose determinant 5(K) / 0. Indeed, if 
substitution of B and ¥ in the second and third of equations ( 33 ) 
is made, the matrix K appears in each term as a common right-hand 
factor, and the equations obtained are therefore consequences of 
equations ( 33 ). 


A final remark concerns some particular solutions of equations (33): 
If D represents a solution of the first of equations ( 33 ) 
and 6 (D) ^ 0, then any scalar multiples of D 


B = ID and ¥ = pD 


with constant factors X and p, p j 0 , satisfy the second and third 
of equations (33)- The proof is immediately obtained by substitution. 

The general solution of system ( 33 ) depends on 12 arbitrary 
constants, but only those solutions which guarantee the preservation 
of property (a) are of interest. By the addition of this new postulate, 
a subclass of the solutions of system ( 33 ) will be obtained. ¥ithout 
trying to characterize the latter completely, a six-parameter 'family 
of solutions of this subclass will be constructed by a specialization 
of the form of the matrices involved. 
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A single method, of assuring the preservation of property (a) is 
to assume ¥ to he of the 0 form 



(34) 


Indeed, the matrix H of the primed equation £' = is supposed 

to have the form ^ 



and the matrix V*lH¥ of the transformed equation 


then 




t s = ¥- 1 HK£ t is 



which shows that property (a) is preserved. 


In order to obtain solutions of system (33) with the matrix W of 
the form of equation (34), an adjustment of the forms of D and B 
must he made hy setting 


(35) 

°) 

It will now he proved that system ( 33 ) has solutions of the form of 
equations (34) and ( 35 ) with arbitrarily preassigned initial values of 

12 12 1 2 

the six available matrix elements 5, 5, 3, 0 , 0 ), and co . 

2 1 2 ' 1 1 2 
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In order to simplify the necessary computations the following 
abbreviations will "be used; 



(35a) 


The 

explicit 


first equation of system ( 33 ) will now he written in an 
form. A simple computation shows that 



( 0 

\ 2 r 1 r 2 
\r) 5 6 


tj 1 6 1 & 2 



the first of equations ( 33 ) is equivalent to the pair of scalar 
equations 


(b X ) s + ^(d 1 ) 2 - t] 1 5 1 6 2 = 0 

( S 2 ) s + t , 1 ^ 2 ) 2 - r,Vs 2 = 0 


(36) 


The integration of the system of differential equations, equa- 
tions ( 36 ), can he reduced to the integration of only one equation with 
one unknown function. As was observed before, the determinant of D, 
given here by 5(D) = -5^6 2 , is constant. This can be verified again 
by multiplying the first of equatipns ( 36 ) by S 2 and the second by 5- 1 - 
and adding the resulting equations. This leads to (s 2 S 2 ) s = 0, which 
shows the constancy of 6 (D) = — o 2 S 2 . If the abbreviation — 6^5 2 = a 
is used, the first equation of system ( 36 ) becomes 
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( sl ) s + , i 2 ( 51 ) 2 + at] 1 = 0 


( 37 ) 


This is a differential equation of first order in 8 1 alone and is of 
the Eiccati type. If a ^ 0 , the second of equations (36) is automati- 
cally satisfied "by S 2 = If a = 0 hut 8 1 j- 0, then 8 2 = 0. 

8 ^ 

All thi s can easily he verified hy the reader. 

1 p 

After equation (37) has been integrated, 3 and 3 are to he 
determined from the second equation of system ( 33 ). A simple computa- 
tion yields 


(s \ 2 - 8 vy 


DHB — HDB = 


(sv _ svy 0 


and the equations sought are therefore 

(p 1 ) s + - &VV = 0 

(p 2 ) s - (®Y - 8 V> 2 - 0 


( 38 ) 


A comparison of system (38) with system (36) shows that in the 
general case where a ^ 0, that is, where both 8 1 and 8 are not 
zero, 3 1 and 8 1 and also 3 2 and 8^ differ only by a constant 
factor, that is. 


3 1 = sV 


2 _2 2 
3 = 8 c 


c 2 , Constant^ 


( 39 ) 
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This is in. agreement; with a result obtained before, which showed that 
in general, if 8(D) ^ 0 , B is obtained from D by multiplying the 
latter by a constant matrix from the right. Indeed, system" (39) can be 
expressed in the form 


B = D 


In the degenerate case where a = 0 and, for example, 8 1 
but 8 2 } 0, equations (39) have to be replaced by 


= 0 


0 1 = — 
82 


p 2 = S 2 c 2 


(°h 


c^ constant 


( 39 a) 




1 2 

If 5 and 5 are both zero, equations (39) are to be replaced 


r 1 1 
0 = c 


P 2 -e 2 


( cl * 


c constant 


( 39 ^) 


as can be easily verified. The essential result, is that, once the 8»s 
are known, the 0‘s can be written down immediately without any inte- 
gration. 
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The differential equations in op- and a> 2 will now he derive d. 
The terms occurring in the third of equations (33) aj ' e 



The differential equations are therefore given by 

(a) 1 ) s - Tj 1 p 2 - T) 1 5 2 (XI 1 + p 1 !] 2 + 5 1 Ti 2 ai 1 = 0 

(cu 2 ) s - T^p 1 - rfb 2 ™ 2 + p 2 ^ 1 + S 2 ^ 2 = 0 


> 


( 40 ) 
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As before, the system of equations (40) can he reduced to one equation 
of the Riccati type in r or r. As was found before, the determi- 
nant S (W) = oAo 2 is constant. This can be verified immediately by 
multiplying the first of equations (40) by o 2 and the second by oi 
and adding * The result is that m 2 + co 2 oo 1 = 0 or that (ar L to 2 ) s = 0 
which implies that oAu 2 = b is a constant. The constant b is 
different from zero because it was always assumed that W is non- 

singular. Replacing a) 2 in the first of equations (40) by — yields 

cc 1 

1 

Ws + ^ ( m1 ) 2 + ( 6l n 2 - S 2 , - T& = 0 (41) 

« 

The second of equations (40) is an automatic consequence of equa- 
tion (4l). This may be verified by replacing cs 2 by to 2 = — in the 

■ aP- 

second equation of system (40). Equation (4l) is of the Eiccati type 
as was asserted before. 


3- Equations Transformable into the Cauchy-Riemann 
Equations or into Those Corresponding to the 
Wave or Tricomi" Equation 

The essential objective of this investigation is to find equations 
of state, or the corresponding density-speed relations, which lead to 
equations transformable into a well-known canonical form. In view of 
this, it is of particular interest to study those equations which can 
be transformed into the Cauchy— Eiemann equations or into those corre- 
sponding to the wave or Tricomi equation by transformations discussed 
in the preceding section. Corresponding to each of these three 
problems, the matrix H of equation (27) becomes, respectively. 



In each case it is required to solve equations (36), (38), and (40), 
where the values of and q 2 from the chosen canonical form hav e 

been substituted. Since all of equations (36), (38), and (40) are of 
the first order, their general solution depends on six arbitrary 
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constants or, in otter words, there exists a six— parameter family of 
transformations leading to the given canonical form. It would he 
incorrect to conclude, however, that the equations transformable into- 
the given canonical form also constitute a six— parameter family. 

Indeed it was already observed that if B- and W are multiplied by 
the same constant matrix from the right, the resulting matrices together 
with the unaltered I) again represent a solution of equations (3 6), 
(38), and (40). In particular this matrix factor may be chosen as a 
constant scalar multiple cl of the unit matrix I. In this case, 
however, (c¥)~ ^-H(cW) = so that the corresponding differential 

equations are identical. This shows that, in general, a five— parameter 
family of equations can be transformed into a given canonical form. 

Consider first the case where the given canonical form is that of 
the Cauchy— Biemann equations. The values of q-*- and tj 2 are 
then q 1 = 1 and t f- = —1. Equations (37) and (4l), to which equa- 
tions (36) and (40) were previously reduced, now become 

5 1 - (s 1 ) 2 + a = 0 (43) 

s ' 


CD^ — ~ + 5 2 )tD^- — 


p 2 = 0 


(44) 


If the constant a is chosen different from zero, the substitutions 

into equation (44) of S 2 = - 4j-, P 1 = S^ 1 , and p 2 = 6 2 c 2 are to 

be made with c^- and c 2 as arbitrary constants. In the degenerate 
case where a = 0, the substitutions to be performed are given by 
equations (39a) or (39b) of the preceding section. 

Equation (43) can be solved in terms of elementary functions, but 
equation (44), it seems, cannot be solved in general in an elementary 
way. In order to obtain solutions expressible in terms of elementary 
functions, the restriction that 5^ = S 2 = 0 will now be made. 


According to equation (39h), P- 1 - 
equation (44) reduces to 


and P^ are then constants, and 


x 2 -o 


4 - x x P) a 
( xl - £ * * 2 - p2 ) 


( 45 ) 


where X-*- and X 2 are arbitrary constants. 
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Equation (45) can be solved, by elementary methods, and the reader 
will easily verify the following statements : 

1 2 

(a) If X and X have different signs, the solutions of equa- 
tion (If 5) have the form 


03 


— A coth a(s — s G ) or 03"*" = A tanh a(s — s Q ) 


(46) 


with arbitrary constants A, a, and s Q . 

(b) If X 1 and have the same- sign, the solutions are of 

the form 


cx?~ = A tan a(s — s Q ) 


(47) 


with arbitrary constants A, a, and s Q . 

(c) If X or X is zero, the solutions have the form 


co 1 = 1 


as + p 


or 03 = as + p 


(48) 


with constant values a and p. 

The computation of the differential equations C s = 1 HW^ t which 
are transformable into the Cauchy-Ri emann differential equations follows 
immediately. A simple computation shows that 



and the equations therefore always have the form 


^s = “"It 
= _ ^t 


o£ 

0>1 


03 


> 


(49) 
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Since = b is a constant, the function cd(s) in system (49) 

takes on one of the following forms: 

CD = C tanh 2 a(s — s Q ) or . cd = C coth 2 a(s — s 0 ) j 


cu = G 


p 

tan cc(s 


(50) 


cd = ±(as + 3) 2 or cd = ± — — 

(as + 3) 2 

where C in the first two of equations (50) is also an arbitrary 
constant. Each formula of system (50) gives a three— or two- 
parameter family of differential equations of the form of equa- 
tions (49) which can he transformed into the Cauchy— Eiemannn 
differential equations. 

Similar considerations lead to differential equations which can he 
transformed into those connected with the wave equation 
(where ql = q 2 = 1). Since these equations can he computed in a 
manner completely analogous to the preceding computation, only the 
final result will he stated: 

All equations of the form 


is = ®nt 
_ 1. 

^S oo^t 


(51) 


with an cd again represented hy one of the functions given hy equa- 
tions (50), can he transformed into equations connected with the wave 
equation. 

The study of equations transformable into a form corresponding to 
the Tricomi equation (where q 1 = 1, q 2 = -s) is of particular interest 
since it can he utilized in the investigation of transonic flows. As 
in the previous discussion, only the case 5^ = S 2 = 0 will he consid- 
ered in more detail. Equation (4l) here takes on the following form: 

CDg — 1 1 s(cd 1 ) 2 \ 2 = 0 

(x* - £ . x* - 0?) 


( 5 2 ) 
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where X 1 and X 2 are arbitrary constants. It seems impossible to 
express the solutions of equation (52) in terms of elementary functions 
except when X, x = 0 or X 2 = 0. 


If 1“^ — the solutions of equation (5 2 ) are given by 


as = as + p 


(53) 


where a and p are arbitrary constants. 

2 

If X = 0 the solutions are given by 


CD 


p 

as + p 


(5*0 


where a and P are arbitrary constants. 

The differential equation £ g = which is transformed into 

equation (27) can here again be computed immediately, and it is seen 
that 


1T 1 HW = 


4 

V° 


0 V ° 



fo 

vs 




2 ' 

CD 1 


7 


The equations therefore have the form 


| s = CDTl t 


1 . m sh 


(*-S) 


(55) 


where cd may take on one of the forms 



30 


MCA TI'T 2065 


03 = ± 


03 


(as + p)2 


= ±(as^ + p)^ 


(56) 


with, arbitrary constants a and p, according as o>*- is given by 
expression ( 53 ) or ( 5 * 0 - 


Equations of the form of equations ( 55 ) may be brought into an, 
at times, more suitable form by a transformation of the independent 
variables given by 


s* = s*(s) 

► 

t* = t 


(57) 


This transformation will, if suitably chosen, simplify the first of 
equations (55) to the form The transformed equations are 


. ds 


_ _ -s ds „ 

a) dS* &t * 


(58) 


The desired simplification is then achieved if 


(R c* 

03 = 1, that is, if 

ds* * 


s* 


J 03(s) ds + Constant 


(59) 


Equations ( 58 ) then take on the form 


where 


^s* ^t* 

* 

= -<D*(s*)| t * 


( 60 ) 


03* ( S*) = 


s 


( 61 ) 
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If as is given by the first of equations ( 5 6 ), the resulting transfor- 
mation is 


s* = ± + Constant- (62) 

3 (as + 0 ) ’ 

A further simple computation, which may he left to the reader, yields 


05 *(s*) = ± P 6 (s* ~ so*) 

[l + ap(s* - s 0 *)]5 


( 63 ) 


where s Q * is an arbitrary constant. If 05 is given by the second of 
equations ( 5 6 ), the form obtained for o>*(s*) becomes more complicated. 
Since the latter form is not used in what follows, it will not be 
written down. 


4. Density-Speed Delations Leading to Systems of Differential 
Equations Transformable into the Canonical Forms of 


the Preceding Section 

If the linearized equation of state p = -- + b is used as an 

P 

approximation to the actual equation of state, then the potential 
function f and the stream function satisfy the Cauchy-Riemann 
differential equations in the holograph plane, after suitable inde- 
pendent variables have been chosen. The results of the preceding 
sections supply a method for constructing' more general equations of 
state, or the corresponding relations between density p and speed q, 
where a transition to a simple canonical form is made possible. This 
is achieved by combining the transformations discussed in the preceding 
sections with suitable simple coordinate transformations in the holo- 
graph plane. 


If the polar coordinates q and Q are used in the holograph 
plane, the equations in 0 and ijr are 

1 


~ 


PI 


9 




(6k) 



32 


HACA TTT 2065 


The local Mach number M and the density p are functions of q 
alone and, according to Bernoulli’s law, 


M 


1 dq 

P cLq. 


(65) 


Equations (64) can be brought into a more suitable form by a 
coordinate transformation of the form 


s = s(q) 


*fc = — Q 


j 


( 66 ) 


It may be noted here that the change in the sign of 9 is introduced 
only to eliminate the minus sign from the first of equations (64). 


The following considerations will first be restricted to the sub- 
sonic region (M <1), where equations (64) are of elliptic type. It is 
well known that equations (64) can then be brought into the symmetric 
form of equations (49) by a suitable transformation of the form of 
equations (66). Indeed after such a transformation, equations (64) 
become 




1 - M 2 do, 


pq ds 
p dq s 


(67) 


Equations ( 67 ) are made symmetric by the choice of a transformation, 
equations (66), which has the property that 


1 — dq _ q ds 

pq ds p dq 


( 68 ) 


This leads to the equation 



(69) 
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or 


= \k ~ M 2 = ]£ + oM 


dq. 


( 70 ) 


which, implies that 


s = 


1 - M £ 


■dq.+ Constant 


( 71 ) 


The function c»(s) which appears in the now symmetrized form of equa- 
tions (67) 


= ^(s)i{r + 


♦b 


( 72 ) 


is given by 


cn(s) 


_ Vl - M 2 


£ 

P 


( 73 ) 


where the right-hand term is to be thought of as expressed in terms 
of Sj as given by equation (71). 

A combination of these results with those of the preceding section 
allows the construction of such density-speed relations which lead to 
equations transformable into the Cauchy-Biemann equations. Each formula 
in system (50) represents a three-parameter family of functions to(s) 
for which equations ( k- 9 ) can be transformed into the Cauchy-Biemann 
equations. If the function cd(s) of equation (73) is now identified 
with a function of one of the families, equations (50), the equation 
thus obtained together with equation (70) represents a system of 
ordinary differential equations of first order in p and s as 
functions of q. By solving this system, with unrestricted values of 
the involved three parameters, a five-parameter family of pairs of 
functions p = p(q) and s = s(q) is obtained. The functions p = p(q) 
are then the ones sought. Since the addition of a constant to s(q) 

(with p(q) left unchanged) leads again to a solution of the differ- 
ential system, a four-parameter family of speed-density relations, 
which leads to equations transformable into the Cauchy-Biemann equa- 
tions, is thus obtained. 
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A constant function is included as a special case of the 
functions co(s) given "by system (50). It is to he expected that the 
functions co(s) = Constant will lead to the density-speed relation 
which corresponds to the linearized pressure— volume relation. This 
will now he verified. 

If, instead of q. 


k = log <1 ( 74 ) 

is int roduced as the independent variable and the constant i is 
called c, equation (73) can he written in the form 


£p _ 

dK 



( 75 ) 


The integration of equation (75) hy the method of separation of vari- 
ables leads to the solution 


P = 


he 

\j<L 2 + 'b 2 


( 76 ) 


where h is a constant, as the reader will easily verify. 

A substitution of function (76) into Bernoulli’s equation finally 
yields 



+ Constant 


( 77 ) 


which is the linearized pressure-volume relation. 


For a general function u>(s) the integration of the system 
consisting of equations (70) and (73) can, hy the elimination of p, 
he reduced to an integration of a differential equation of second order 
in s(q) alone. Indeed, equation (73) shows that 


ds 

p =^- 

cd(s) 


( 78 ) 


which, when substituted into equation (70), leads to 


d^s 

dK 2 




dcu 


- 


co( s ) ds \dK/ 


ds 

dK 


( 79 ) 
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after a simple computation. If any solution of equation (79) is substi- 
tuted into equation (78 ) , a desired function p = p(q) is obtained. 

Equation (79) 'becomes essentially simpler if the roles of the inde- 
pendent variable k and the dependent variable s are interc hang ed 
If the abbreviation 


T 


dK 

ds 


(80) 


is introduced into equation (79) t a simple computation yields 


_ T 2 _ d log _ 
ds ds 


(81) 


Any solution of the Eiccati type equation, equation (8l), gives a 
possible function k = k(s) when integrated. 


The function co corresponding to the real equation of state has 
a special analytic character, which is useful to know in the discussion 
of approximating functions of a>. The pressure-density relation of gas 
dynamics, upon which <n(s) depends, is of the form p = ap 7 + fc, 
where a, b, and 7 are constants and 7^1. If dimensionless 
variables are introduced such that the density and velocity of sound 
are both 1 at q = 0, the following formulas are obtained: 


and 





(7 = 1) 



( 82 ) 

(83) 

m 


In addition, the function cu is given by 

..i 1 " 


7 + 1 2 V / 2 


7+1 

1 - 7 ~ 1 q g N ) 2 (^'~ 1 ) 


(7 f 1) 


( 85 ) 


) 


2 
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and 


oi = \j 1 — q 2 e 2 


(7 = 1 ) 


These formulas show that all three functions 
analytic functions of q 2 and can he represented 
which converge in the whole subsonic region q 2 < 

the power series for 01 begins with the terms 


p, M 2 , and oi are 
hy power series in q 2 
— ^~l m 1x1 particular. 


. 7 + 1 li 

oi = 1 - — — - q^ + . . 


( 86 ) 


Formula (71) shows further that if a is defined as 


or 


= e 2s 


(87) 


then the latter function can also he represented hy a power series 
in q 2 . if the constant of integration in equation ( 71 ) is suitably 
chosen, the power series for a begins with the terms 



( 88 ) 


From a combination of equations (86) and (88), it follows that oi can 
he represented hy a power series in a starting with the terms 


oi = 1 — 2 1 . t— 1 . o 2 + . . . ( 89 ) 

8 

The question of how far the function oi(s) can he approximated hy 
a function of the family, equations ( 50 ), can now he discussed. If the 
approximation is desired for a range of q starting with q = 0, the 
best possible adjustment at q = 0 is achieved with a func- 
tion o>l = C tanh 2 a(s — s Q ), with the following choice of the constants: 


C = 1 
a = 2 


4e~ ^ s ° = 2L±_i 


8 


(90) 
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This shows that 



( 91 ) 


The speed— density relations which correspond to a) = ay, have now 
to be studied. In order to find these relations,, it is necessary that 
the function be substituted for o> in equation (8l). If a solu- 

tion of the resulting differential equation is integrated with respect 
to a a possible s-k relation is obtained. A desired speed-density 
relation is then obtained by a subsequent substitution into 
equation (78). 

It is convenient here to use the quantity cr = e 2s of equa- 
tion (87) as the independent variable instead of s. Equation (8l) 
in which m has been replaced by then takes on the form * 


2a 


dv 

dtr 


t 2 + 


2L±la2 

2 


-T - 1 


1 - 


V + 1 ) 

K 32 ) 


( 92 ) 


That particular solution t 1 which allows for the best approxi- 
mation to the "real" speed-density relation equation (82) or (83), 
at q = 0 is now to be determined from among the infinitely many solu- 
tions of equation (92). 

It can easily be deduced from equation (88) that the function t 
which corresponds to the real equation of state has a power— series 
expansion in a starting with the terms 


T 



( 93 ) 


It will be shown in the appendix that equation (92) has one and only 
one analytic solution in the neighborhood of a = 0 which satisfies 
the initial conditions 
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t = 1 



for a = 0 > 


(9k) 


where is an arbitrary, constant. 

The particular solution of equation ( 92 ) corresponding to ^ 

therefore gives the best possible adaptation to the real equation of 
state at q = 0. The s — k relation is finally obtained by an inte- 
gration of this t with respect to s. The constant of integration 
has to be adjusted such that f \ = 1 since tliis is the case with 

\dq2/q = o 

the real equation of state as shown by equation ( 88 ). 

The discussed approximation method is not the only possible one. 
One ma y* for instance* introduce the requirement that the approximating 
function cni(s) coincide with the function cu* which iB deduced from 
the real equation of state* at three preassigned points. Numerical 
results of the mentioned approximations will be discussed at the end of 
this section. 

The foregoing investigation has so far been restricted to subsonic 
flows* and will now be extended to the transonic region.- The canonical 
form* equations (t 9 ) * must now be abandoned since it can arise only for 
an elliptic— type system of differential equations. This form will be 
replaced* as is usual for a mixed flow* by a new one for which the 
second equation of system (6j) becomes 

- -* s (95) 

This is achieved by a transformation* equations ( 66 )* which satisfies 
the condition that 


^=2 ( 96 ) 

dq q 


or 



+ Constant 


(97) 


The first equation of system (6j) takes on the form 

0 S = <»*(s>jf t 


(98) 
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with. 


ov*(s) 


1 - M 2 


p 4 cl 


(99) 


In the case where the speed-density relation,, equation (82), is used, 
then co*(s) becomes 


or 


o)*(s) = 




7+1 

i-^r 1 


(7 i 1) 


co*(s) = (l - q 2 )e^ 2 (7 = 1) 


( 100 ) 


( 100 a) 


Remark.— A comparison of equations (100) and ( 85 ) shows that in 

the subsonic region and for equal -values of q the function ay* is 
the square of a>. This relation could also have been obtained by a 
consideration of invariants of the differential equations under coordi- 
nate transformations. 

The following procedure is now parallel to that followed in the 
subsonic region. It is desired to construct speed— density relations 
which lead to equations transformable into those connected with the 
Tricomi equation. This is achieved by identifying the function o>* of 
equation ( 99 ) with a function of formula ( 63 ) in which s* is now 
replaced by s. Equations ( 96 ) and ( 99 ) represent then a system of two 
ordinary differential equations of first order in p and s as func- 
tions of q or k = log q. -Again this system can be reduced to a 
differential equation of second order in s as a function of k by 

substituting p = — from equation ( 96 ) into equation (99). The result 
is cue 

. 0 

dn 2 \dK/ 


( 101 ) 



4o 


MCA TN 2065 


'Tills again is simplified by choosing s as the independent variable. 

cLic 

If the substitution of r = — is performed as before, equation (101) 

ds 


reduces to 


— = t ^ s ) ( 102 ) 

ds 

/ 

Should t = ^ be used as the dependent variable, then equation (102) 
would reduce to 


dr 

ds 



) 

(102a) 


The question of choosing the most suitable approximation func- 
tion oj*, call it 03 *-l, from among the functions of formula ( 63 ), will 
now be discussed. An approximation is desired in the neighborhood of 

the sonic speed q„ = 1 / . Without loss of generality it- can be 

V7 + 1 

assumed that the value of s corresponding to this sonic speed is zero. 
The constant s Q in formula ( 63 ) must then also be set equal to zero. 
The best approximation of the function cd*(s) of formula (100) 
by a 3 *!(s) will be achieved if the remaining two constants a and 0 
are so chosen that 


/ (kl3*\ _ /dU3*-A 

W ) s=0 \ ds y s=0 



(103) 


By the use of formulas (100), (97), and (82), a simple computation 
yields 




= -2(2 7 


+ 5 ) 


7+3 



( 104 ) 
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A comparison with formula (63) written in the form 


05* 


1 


as 

(1 - bs) 5 


shows that a and b are 



( 105 ) 


(106) 


After the computation of C5*q ( s ) has been completed, equation (102) 
must be solved. The initial condition to be satisfied is obtained from 
the condition that the relation between s and q and the relation 
obtained from the real equation of state must agree at the sonic 
speed q = q s up to the second order. 


Equation (97) shows that 


T 


1 1 

p (1-^)^ 


(107) 


and the initial condition for equation (102) is obtained by a substi- 
tution of q 2 = q 2 = — - — . This yields the condition that 
s 7 + 1 


1 


T = ( 2 L lp ) 7 1 for 5=0 


(108) 


By an integration of the corresponding solution of equation (102) 
with respect to s, the function k = k(s) is obtained. The constant 
of integration must again be adjusted to the real equation of state. 
This requires that 



7 + 1 


for s = 0 


(109) 
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The previous considerations contain constructions of approximating 
speed-density relations in tlie subsonic and transonic regions leading 
to simplified equations of motion. A similar procedure would lead to 
approximation formulas in a pure supersonic region. (Transformations 
which lead to the canonical form connected with the ware equation would 
then have to be used. Since the discussion in such a case would be 
absolutely parallel to the one for the subsonic region, it is omitted 
here. 


In order to have a picture of the approximations which can be 
achieved, numerical computations have been made in the subsonic and 
transonic cases. The results are given in tables 1 and 2 and illus- 
trated in figures 1 and 2, respectively. 

Table 1 contains the tabulation of four functions cd, ci^, oog, 
and ou. The first one is the function given by formula (85) 
with y = 1 . 405 . The variable 


°1 


= 1(7 + 



( 110 ) 


is chosen as the independent variable instead of 0 . The variable 0 ^ 
is normalized by the chosen constant factor so that it ranges in the 
total subsonic region from 0 to 1 . The values of CD are taken from 
reference 4 . All the functions cdj_, ci£, and cdo are approximating 
functions of cd of the form of the first of equations (50) and all can 
be written in the form- 


cA-o^y 

V 1 + “«i7 


(in) 


The function cd-^ is the approximating function given by formula (91), 
and it gives the best possible approximation to cd at q = 0 . The 
functions cn£ and cdo are interpolation functions of the form of 
equation ( 111 ) determined by the conditions 


cd = cDg for 0 1 =O(M=O), 02= 0.63(M= 0.4702), and 0 1 = 0 . 80 (M= 0 . 6344 ) 

(112a) 


cd=cd 3 for 0 1 =O(M = O), a ± = 0.63(14= 0.4702), and 0 = 0.96(M= O.8698) 

(112b) 
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In the range, for example, from crp = 0 to a 1 = 0.80, which corre- 
sponds to a range of M from 0 to O. 63 H, the deviation of ov> from ay 
is never greater than one unit of the third dec imal place/ with the 
exception of the last entry where the deviation becomes three units of 
the third decimal place. 

Table 2 contains the tabulation of the function ay* given by 
formula ( 100 ) which occurs in the discussion of the transonic region and 
the approximation function oy*j_ given by formula ( 105 ), for the 
value y = 1 . 5 . The approximation of ay* usually used is that given 
by the linear function oy* 2 (s) giving the tangent of the ay*-curve at 
the zero point which corresponds to the sonic speed. For comparison, 
all three functions ay*, ay*^, and oy * 2 are represented in figure 2 . 

It is seen that down to s = —0.22, which corresponds to a Mach number 
of about M = 0.70, oy*-^ gives a far better approximation to ay* 

' than oy* 2 . Also in the supersonic region up to, for example, s = 0.11, 
which corresponds to a Mach number M = I. 29 , the function ay* is 
strongly superior to ay* 2 . ^ 


CONCLUDING REMARKS 


The present investigation contains the study of a class of trans- 
formations of the Baecklund type which transform systems of partial 
differential equations into each other where simple point transformations 
ordinarily fail. . Corresponding solutions of the two 'systems are derived 
from each other by the solving of a system of ordinary differential 
equations, a process, which, both from a theoretical and practical point 
of view, is essentially simpler than solving a system of partial differ- 
ential equations. An application of this study to gas dynamics leads 
to a five— parameter family of speed-density relations for which the 
corresponding equations for the stream and potential functions can be 
transformed into the Cauchy-Riemann equations or into those connected 
with the wave or Tricomi equation. 

The methods developed allow for a far-reaching extension explained 
as follows. 

The class of transformations studied does not form a group; that 
is, a transformation obtained by a composition of two such transfor- 
mations is in general not contained in the original class of transfor- 
mations. This represents an advantage, since the process of composi- 
tion can be repeated arbitrarily often, and the class of transformations 
can be extended more and more. This iteration process allows the 
construction of a family of speed— density relations, which contain 
arbitrarily many parameters and which lead to equations transformable 
into one of the three canonical forms used. 
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It is of great theoretical and practical interest to determine all 
systems of partial differential equations whose form can he brought 
arbitrarily close to one of three canonical forms by such composite 
transformations . 


Syracuse University 

Syracuse , N. Y . , October 27, 1947 
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APPENDIX 


PROOFS OF THEOREMS 


The theorem used in section 4 will he proved in the following 
generalized form: 

Theorem 1 . Consider the differential equation 

2d ^ = t 2 + k(a) r - l . (Al) 

where the coefficient k(cr) is an analytic function of cr in the 
neighborhood a = 0 and has a zero at cr = 0 of at least the 
second order. The power— series expansion of k(cx) is therefore 
of the form 

k(cr) = a^ 2 + a^cr^ + . . . (A2) 

If k(cr) satisfies these conditions,, then there exists a solution 
of equation (Al), analytic about a = 0, which satisfies the 
initial conditions 



(A3) 


where c-^ is arbitrary. 

Proof .— The required solution must be expansible in a power series 
of the form 


T 


= 1 + c^cr + 


) a 2 


(Ah) 



h6 
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By substituting series (Ah) into equation (Al) and comparing coeffi- 
cients of equal powers on both, sides of equation (Al), the relation 


2(n - l)c n 


n— 1 


n— 1 


c i c n— i + 2_ a i c n— i + «n 
i=l i=l 


(A5) 


(n — 2, 3j • • • ) 


is obtained, where a^ = 0. 

Equation (A 5 ) represents a recurrence formula which allows for a 
successive confutation of all coefficients Cj_ in series (Ah) . It will 
be proved now that series (Ah) has a positive radius of convergence. 
Series (A 2 ) has a majorant of the form 

k(tf) = = a^a 2 + a^ff 3 + . . . (A6) 

having a radius of convergence which can be assumed arbitrarily close 
to that of series (A 2 ). 

The values of the coefficients in series (A 6) are 

= Aa 1 “ 2 (i = 2 , 3, . . .), = 0 (A 7 ) 

If in formula (A 5 ) the a^’s are replaced by the a^*s "the 

factor n — 1 on the left side is dropped, a modified recurrence 
formula 


n— 1 n— 1 



2 °n = 2Z °i°n-i + 21 ^i^n-i + % 
i=l i=l 

( n - 2, 3, . . . ) 

(a8) 


is obtained, which again determines c n uniquely once c^ has been 
chosen. 

If a value for c^ is chosen such that 




(A 9) 
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then a comparison of formulas 


(A5) and (a8) 


shows immediately that 



(n = 1, 2, . . . ) 


(A10) 


It is, however, easy to show that the series 


c-j_cr + C2(j2 + 


(All) 


has a positive radius of convergence. Indeed just as the c n *s, given "by 
recurrence formula (A5), are the coefficients of a solution to the 
differential equation, jequati on (Al), equation (a 8) expresses the 
requirement that the c n 's are the coefficients of a generating 
function 


00 

T(<r) = YZ ^ 


n=l 


(A12) 


which satisfies the algebraic relation 


or 


2t = T 2 + - ^ — (1 + t) 

1 — CtO 2 




(A13) 


(AHO 


and which vanishes for a = 0. This determines 7 to he 


A0 2 J 

Ti a ° 2 

2 2 
AcT 

2(1 - Oc 2 ) y 

2(1 - aa 2 )_ 

1 — acr 2 


(A15) 


But function (A15) can he expanded into a power series in a with a 
positive radius of convergence. This completes the proof of theorem 1. 
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It is of interest to know how far the solution t(ct) of equa- 
tion (Al) satisfying initial conditions (A3) is represented "by the 
power series , series (A4). A definite answer can he given under the 
restrictive conditions 


tp 
11 V 
0 

(n — 2, 3j • • • ) 

(Al6) 

c 1 ^ 0 


(Al?) 


It will he shown later that these conditions are fulfilled when equa- 
tion (8l) contains an to corresponding to the real equation of state, 
provided that 7 = 1, or when it contains an <x> = coj_, a function which 
was used as an approximation to the former and given hy equation (91). 

The following theorem is given in answer to the foregoing question; 

Theorem 2 . If the restrictive conditions (Al6) and , ( A17 ) are 
satisfied and one knows that the solution r exists in a range 
0 < a < a QJ then the radius of convergence of series (a 4) is at 

least equal to cr 0 . 

Proof.— Now recurrence formula (A5) shows immediately that 

c n > 0 (n = 1, 2, 3, . . .) (Al8) 

An immediate consequence of inequalities (Al8) is that if a 1 represents 
the radius of convergence of series (A4), all derivatives of t(ct) in 
the interval 0 ^ cr < oq 8X0 non-negative. An analytic continuation 
hy power series shows further that this is true not only for 0 < a < 
hut in the whole interval 0 S a < a Q . 

A theorem of S. Bernstein (reference 12) is now used to complete 
the proof. Bernstein’s theorem states; 

If a function t(c) has derivatives of arbitrarily high 
order and al 1 are non-negative in the interval 0 % a < a Q (the 
function is then said to he absolutely monotonic in the interval) 
then it is analytic there and can he represented hy a power series 
in a whose radius of convergence is at least equal to a Q . 
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In the special cases considered in section 4, the value of cj_ 
was equal to Condition (A17) is therefore satisfied. But condi- 
tion (Al 6 ) is also fulfilled. This is expressed in 

Theorem 3 . . Condition (Al 6 ) is satisfied "both 
for k(s) = — — log a), where 03 is represented hy formula ( 85 ) 

with 7^1 and for ^(s) =-|- log 03j_ with 0 ^ given hy 
formula ( 91 )- 

Proof .— That the theorem is true for kj_ is immediately seen from 

the expression for in equation (92). In order to verify it 

for k(s), it must he shown that the function k and all its deriva- 
tives with respect to a are non-negative. 


Consider first k itself: 


k =• 


d log 03 _ d log 03 dq^ 


ds 


dq^ 


ds 


(A19) 


According to equations ( 85 ) and (71), 


d log 03 _ 7 + 1 
= 1 


(x _ 1^(1 - ^ 2 ) 


2 /I - 

^ = 2<lV 2 

48 ' 1 - XTi 2 


(.A20) 


(421) 


Functions (A20) and (A2l) have power-series expansions in q 2 with non- 
negative coefficients. This is immediately clear for function (A20). 

In order to prove it for function (A2l), consider 


log | 


8 

OJ 

p 

H 

1 

N 

1 

H 

(7 + iy ( 


2 _in 

7 - i\ n 

1 7 ± lg2 2 n 

V 2 ) \ 

2 7 . 


, 2 n 


(A22) 


2 
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All the coefficients of series (A 22 ) are certainly non-negative 
since 7 ^ l s and therefore those of the expansion of function (A21) 
and those of equation (A 19 ) are also non-negative . Thus the fraction k 
is non-negative. 

It vill now he shown that 


dk _ dk dq 2 1 
da dq 2 ds 2a ' 


(A 23 ) 


is also non-negative. The first two factors of the right-hand member 
of equation (A23) have power-series expansions in q2 with non-negative 

coefficients. It mil now he shorn that — , which is obtained as a 
o . . a 3 

function of q from equation (71) j has an expansion of the form 



(A 24 ) 


all of whose coefficients 7 n are 
prove that all the coefficients of 


non-negative . 
the expansion 


It 

of 


is sufficient to 



are 


non— positive with the exception of 


l_2U-kq2 


= 1 - 


1 - 2^2 

2 


n=l 


the first. Indeed, 

(LjrT 1 1 2 ” - 1 - ti 


where all the coefficients of the expansion for U are non-negative. 
Now 



CO 


n=0 


(-D 



2 

evidently involves only non— positive coefficients in q with the 
exception of the first term 1 . By a multiplication of the series 
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In the special cases considered in section A, the value of cj_ 
was equal to Condition (A1J) is therefore satisfied. But condi- 
tion (Al 6 ) is also fulfilled. This is expressed in 

Theorem 3 . . Condition (Al 6 ) is satisfied "both 
for k(s) = — — log <u, where o> is represented by formula ( 85 ) 

.with 7^1 and for k-^s) =-^- log with given hy 

formula ( 91 )- 

Proof .— That the theorem is true for k^ 

the expression for k-j. in equation ( 92 ). In 
for k(s), it must "be shown that the function 
tives with respect to a are non-negative . 

Consider first k itself: 


is immediately seen from 

order to verify it 
k and all its deriva— 


k d log ca d log pi dq 2 
ds dq 2 ds 

According to equations ( 85 ) and (71), 


d log 03 _ 7 + 1 
dq 2 

d-Q. 2 _ 
ds 




(A19) 


(A20) 


(A21) 


Functions (A20) and (A2l) have power-series expansions in q 2 -with non- 
negative coefficients. This is immediately clear for function (A20) . 

In order to prove it for function (A2l), consider 


log, 


1 1 

00 

r 


2 - 

- 1 y~ 1 

( 7 + lY 5 * / 

7 - i\ n 

ll 

2 5=1 n 

\ 2 J V 

2 / . 


, 2 n 


(A22) 


2 
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All the coefficients of series (A22) are certainly non-negative 
since 7 ^ 1, and therefore those of the expansion of function (.A21) 
and those of equation (A19) are also non-negative. Thus the function k 
is non-negative. 

It will now he shown that 


dk _ dk dq 2 1 
dcr dq 2 ds 2 a 


(A23) 


is also non-negative. The first two factors of the right-hand member 
of equation (A23) have power— series expansions in q^ with non-negative 

coefficients. It mil now he shown that which is obtained as a 

function of q 2 from equation ( 71 ) , has an expansion of the form 


- = -^(l + 7 X 1 2 + 7 2 1 4 + • • -) 


(A24) 


a 1 1 of whose coefficients 7 are non-negative. It is sufficient to 

fl - — - q 2 

prove that all the coefficients of the expansion of | j are 

1 - 

2 

non-positive with the exception of the first, indeed. 


1 -£+_iqa 


l-£— iq 2 


-i-f: i-v 

n=l \ ^ 7 


where all the coefficients of the expansion for U are non-negative. 
Now 

2 


7 — 1 2 
1 — — q^ 


= ]f l - u = XI (-in 2 

n =0 ' 


evidently involves only non— positive coefficients in q with the 
exception of the first term 1. By a multiplication of the series 


MCA TIT 2065 


51 


expansions of the factors in the right-hand member of equations- (A23) ^ 

the result is obtained that ^ has an expansion in q2 with non- 
negative coefficients and is therefore non-negative. 

tS 

By a continuation of these considerations in the same manner it 
can evidently be proved that all derivatives of k with respect to a 
have power-series expansions in q 2 with non-negative coefficients and 
are therefore non-negative. 



52 


MCA TN 2065 


DEFERENCES 


1. Chaplygin, S. A. : On Gas Jets. Scientific Memoirs, Moscow Univ- , 

Math. Phys. Sec., vol. 21, 1902, pp. 1—121. (Eng. translation 
published "by Brown Univ. , 1944; also NACA TM IO 63 , 1944* ) 

2. Ton Ear man, Th. : Compressibility Effects in Aerodynamics. Jour. 

Aero. Sci., vol. 8 , no. 9 , July 1941, pp. 337— 3 56. 

3 . Tsien, Hsue-Shen: Two-Dimensional Subsonic Flow of Compressible 

Fluids. Jour. Aero. Sci., vol. 6 , no. 10, Aug. 1939, PP* 399-407-. 

4. Bers, Lipman: On a* Method of Constructing Two-Dimensional Subsonic 

Compressible Flows around Closed Profiles. NACA TN 969 , 1945* 

5 . Bers, Lipman: On the Circulatory Subsonic Flow of a Compressible 

Fluid past a Circular Cylinder. NACA TN 970, 1945* 

6 . Bers, Lipman: Telocity Distribution on Ving Sections of Arbitrary 

Shape in Compressible Potential Flow. I — Symmetric Flows 
• Obeying the Simplified Density-Speed Eelation. NACA TN 1006, 1946. 

7. Gelbart, Abe: On Subsonic Compressible Flows by a Method of Corre- 

spondence; I — Methods for . Obtaining Subsonic Circulatory 
Compressible Flows about Two-Dimensional Bodies. NACA TN 1170, 
1947. 

8 . Bartnoff, Shepard, and Gelbart, Abe: On Subsonic Compressible Flows 

by a Method of Correspondence. II — Application of Methods to 
Studies of Flow with Circulation about a Circular Cylinder. 

NACA TN 1171, 1947. 

9 . Lin, C. C. : On an Extension of the Ton Karmsm-Ts i en Method to Two- 

Dimensional Subsonic Flows with Circulation around Closed Profiles. 
Quart. Appl. Math., vol. IT, no. 3 , Oct. 1946. 

10. Darboux, Gaston: Lecons sur la th^orie gdnerale des surfaces. 

T. III. Gauthier-Tillars (Paris), 1894, pp. 438-444. 

11. Weyl, Hermann: Mathematische Analyse des Eaunproblems . J. Springer 

(Berlin), 1923, PP* 66-68. 

12. Widder, David Ternon; The Laplace Transform. Princeton Math. Ser., 

vol. 6 , Princeton Univ. Press, 1941, p. 146. 



NACA TI'I 2065* 


53 


TABLE 1.— VALUES OP to, a>j, <n>, AUD tOj 


t 


to given by equation (85), where — = e s , relationship "between s anfl 

1.45 


q given 


(l - off P\ 

"by equations (71) and (84), and 7 = 1.405. <n> and ctb obtained from [ i-] 

\1 + °°i / 

which is an approximation formula of the first of equations (50). Constants a 
and fi evaluated for = O.63 (Mach number, 0.4702) and = 0.80 (Mach 
number, 0.6344) in teg and for tr^ = 0.63 (Mach number, 0.4702) and cr^ = O.96 
(Mach number, O.8698) in to^.j 


*1 

CD 


“2 

“3 

0.02 

1.0000 

1.0000 

1.0000 

1. 0000 

.16 

1.0000 




.18 

•9999 




.20 

■ 9999 




.22 

.9998 

1.0000 



.24 

■ 999 8 

.9996 

.9999 


.26 

•9997 

.9996 

.9999 


.28 

.9996 

.9996 

.9998 


• 30 

•9994 

.9996 

.9998 

■ 9999 

.32 " 

.9992 

.9992 

.9996 

-9999 

-34 

•9990 

.9992 

.9994 

.9998 

•36 

.9987 

.9988 

.9992 

-9996 

.38 

.9984 

.9988 

.9990 

-9996 

.40 

.9980 

.9988 

.9986 

-9994 

.42 

•9975 

.9980 

.9982 

-9990 

.44 

.9968 

.9976 

. -9978 

.9988 

.46 

.9961 

.9968 

.9972 

.9982 

.48 

-9953 

.9984 

.9964 

-9976 

.50 

.9944 

.9956 

-9954 

.9968 

•52 

•9933 

• 9948 

.9942 

-9956 

• 54 

- .9920 

.9940 

.9930 

.9944 

.56 

.9905 

.9932 

■9914 

-9926 

.58 

.9888 

.9924 

.9894 

-9906 

.60 

.9868 

.9912 

.9872 

.9878 ' 

.62 

.9847 

.9900 

.9849 

-9947 

.64 

.9818 

.9884 

.9819 

.9809 

.66 

.9791 

.9872 

-9785 

.9761 

.68 

-9749 

.9853 

.9746 

-9706 

.70 

-9709 

.9837 

.9702 

.9839 

•72 

.9660 

.9817 

.9651 

-9559 

.74 

.9602 

.9799 

.9596 

-9465 

• 76 

-9535 

.9775 

• 9532 

-9357 

.78 

-9455 

.9752 

.9460 

.9229 

.80 

-9350 

.9724 

.9380 

.9082 

.82 

.9238 . 

.9696 

.9291 

-8913 

.84 

-9113 

.9667 

• 9193 

.8720 

.86 

.8948 

.9631 

.9086 

.8503 

.88 

.8740 

.9596 

.8966 

.8257 

.90 

.8492 

.9561 , 

.8836 

.7987 

.92 

•8153 

.9520 

.8694 

.7681 

• 94 

.7720 

.9479 

.8540 

-7350 

.96 

.7037 

.9436 

.8372 

.6984 

.98 

.5891 

.9388 

.8192 

.6595 

1.00 

.0000 

.9339 

.7998 

.6173 
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TABIE 2.- VALDES OF CD* GIVEN BY FCEMOLA (100) 


AND cd*^ GIVEN BY FQBMDLA ( 105) 


|ln formula (100), s = J — dq; 

formula (105), a and "b de 
equations (106) for 7 = 1-5- 

y = 1.5. In 
jermined "by 

s 

CD* 

cd*2 

-0.30 

0.8832 

0.5821 

-.28 

.8687 

• 5955 

-.26 

.8520 

.6072 

-.24 

.8329 

.6165 

-.22 

.8108 

.6227 

-.20 

.7849 

.6250 

-.18 

.7548 

.6223 

-.16 

• 7195 

. 6132 

-.14 

.6779 

.5962 

-.12 

.6283 

.5690 

-.10 

.5690 

.5292 

-.09 

• 5349 

.5037 

-.08 

• 4973 

.4738 

-.07 

.4559 

.4389 

-.06 

.4101 

.3986 

-.05 

• 3597 

.3522 

-.04 

.3028 

.2989 

-•03 

.2394 

.2380 

-.02 

.1686 

.1686 

-.01 

.0895 

.0897 

.00 

.0000 

.0000 

.01 

-.1021 

-.1016 

.02 

-.2183 

-.2164 

•03 

-.3503 

-3464 

.04- 

-.5034 

-4930 

.05 

-.6789 

-.6585 

.06 

-.8853 

-.8449 

.07 

-1.1270 

-3.0552 

.08 

-1.4127 

-1.2921 

.09 

-1*7536 

-1.5588 

.10 

-2.1659 

' -1.8594 

.12 

-3.2931 

-2.5792 

.14 

-5.0852 

-3-4935 

.16 

-8.1820 

—4 . 6564 

.18 

-14.2794 

-6.1398 

.20 

-29.0994 

-8.0379 








56 


HACA TN 2065 



NACA-Langley - 4-13-50 - 1000 


Figure 2.- Results of computations in transonic case. 



